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ON THE DESIGN-CONSISTENCY PROPERTY OF 
HIERARCHICAL BAYES ESTIMATORS IN FINITE POPULATION 

SAMPLING 

By p. Lahiri^ and Kanchan Mukherjee 

University of Maryland, College Park and University of Liverpool 

We obtain a limit of a hierarchical Bayes estimator of a finite pop- 
ulation mean when the sample size is large. The limit is in the sense 
of ordinary calculus, where the sample observations are treated as 
fixed quantities. Our result suggests a simple way to correct the hier- 
archical Bayes estimator to achieve design-consistency, a well-known 
property in the traditional randomization approach to finite popula- 
tion sampling. We also suggest three different measures of uncertainty 
of our proposed estimator. 

1. Introduction. Ericson [4] put forward a subjective Bayesian approach 
in finite population sampling. The subjective Bayes and the more general 
hierarchical Bayes estimators have been found to be eff'ective in combining 
information from a variety of sources in conjunction with the sample survey 
data on the main variable of interest; see [10]. Specific applications include 
small area estimation, estimation from longitudinal surveys, and so on. 

In Section 2, we obtain the mathematical limit of a general version of 
Ericson's subjective Bayes estimator of a finite population mean when the 
sample size is large. By mathematical limit, we mean the limit in the sense of 
ordinary calculus, where observations are treated as fixed real numbers. We 
show that the Bayes estimator converges to a quantity which is free of any 
hyperparameters, but may depend on the sample observations, irrespective 
of the model used to derive the Bayes estimator. This pure mathematical 
result can be used to examine the design-consistency property of the sub- 
jective Bayes estimator. Design-consistency is a desirable property in the 
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randomization approach to finite population sampling. Kott [13] advocated 
the use of a model-based estimator of a finite population mean which is 
also design-consistent. For a formal definition of design-consistency, see [22], 
page 18. We find that the subjective Bayes estimator is, in general, not 
design-consistent. In other words, the subjective Bayes estimator does not 
converge to the true finite population mean as the sample size becomes 
large. Here the convergence is defined with respect to the sampling design, 
observations for all units of the finite population being treated as fixed non- 
random quantities. The result, in turn, suggests a simple correction to the 
Bayes estimator to achieve design-consistency. 

In Section 3 we consider the case where the finite population is divided 
into several strata and a sample is available from each stratum. Model-based 
estimation (both empirical and hierarchical Bayes) of a stratum mean has 
been considered by several authors; see Ghosh and Meeden [9], Ghosh and 
Lahiri [6], Battese, Harter and Fuller [2], Datta and Ghosh [3], Prasad and 
Rao [19], Ghosh and Lahiri [8], Arora, Lahiri and Mukherjee [1] and Jiang, 
Lahiri and Wan [12], among others. Prasad and Rao [20] observed that 
these model-based estimators are typically not design-consistent unless the 
sampling design is self-weighting within each stratum, and they proposed a 
pseudo-EBLUP design-consistent estimator for the normal case; see [23] for 
a pseudo-hierarchical Bayes version of Prasad and Rao [20]. However, it is 
not clear how one can extend their approach to a general nonnormal case, 
since the distribution of a linear combination of observations used in their 
paper is not always analytically tractable for nonnormal cases. In Section 3 
we provide the stratified sampling extension of the result of Section 2 for a 
general hierarchical Bayes estimator of the finite population mean when the 
hyperparameters are assumed to be unknown. Again, our mathematical limit 
result lends itself to a simple correction of the hierarchical Bayes estimator 
to achieve the design-consistency property. 

For small area problems with binary data, Folsom, Shah and Vaish [5] 
proposed a method which appears to work well in their simulation work 
in terms of design-consistency. However, a formal proof is needed to claim 
design-consistency of their estimator. In this connection, we also mention 
the work of Malec, Davis and Cao [16], who considered sample selection 
adjustment to their model-based procedure. See [21] and [11] for further 
discussions of design-consistency in small area estimation. The simple and 
general approach taken in this paper merely involves finding the mathemat- 
ical limit of a hierarchical Bayes estimator. 

Measuring uncertainty of the proposed design-consistent model-based es- 
timator is an important problem. In Section 4 we first present two possible 
ways to measure the uncertainty of the proposed design-consistent model- 
based estimator. We note their merits and disadvantages and propose a third 
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measure that is just a simple average of the two. The results from a sim- 
ulation study are presented in Section 5. In our simulation experiment the 
proposed measures perform reasonably well and the third measure appears 
to be a compromise between the first two. The technical derivations of the 
results of Sections 2 and 3 are relegated to the Appendix. 

2. The Bayes estimator and its limiting behavior. Let Yi be the value 
of the characteristic of interest for the ith unit of a finite population (i = 
1,...,N). We assume that the finite population size N is known. In our 
subjective Bayesian formulation, we assume that Yi, i = 1, . . . , N, are in- 
dependent realizations from a superpopulation belonging to the family of 
densities 

(2.1) f{y\e, ((>) = eM4>"'{yO - m} + piy, <t>)], 

where <j)> 0. When (p is known, the generalized linear models characterized 
by the above class of densities include the exponential family of distribu- 
tions (in particular, Gaussian, Bernoulli and Poisson distributions). Thus, 
it includes continuous data as well as both the categorical and count data. 
When (p is unknown, it includes distributions which are not in the exponen- 
tial family. Note that the function ip is given by 

^{e)=(j)log J exp[(l)~^ye + p{y,(l))]dy 

and the mean of Y is ip'{9); see [18], page 28. Here and throughout the 
paper, a' denotes the derivative of an arbitrary function a. 

We assume that the superpopulation parameter 9 has the prior distribu- 
tion 

(2.2) h{e) = p + u, 

where h, called the link function, is a strictly increasing function of 9, (3 GM 
is a location parameter and u~ N{0,r~^). Similar Bayesian models were 
used by Ghosh and Meeden ([10], page 269). 

A sample of fixed size n is drawn from the finite population. Let {p{s)} 
denote the sampling design. Note that p{s) is the probability of drawing a 
particular sample s of size n from the universe of all possible samples S of size 
n. Thus, p{s) > for all s € 5 and J2sesPi^) ~ ^- Let VTj = P{s B i) be the 
first-order inclusion probability of unit i, that is, the probability of including 
the population unit i in the sample (i = 1, . . . , N). Let ys = {Yi,i G s) be the 
vector of observations in the sample. 

The standard Bayesian approach to finite population sampling recognizes 
the importance of a good sampling design for selecting the sample, but once 
the sample is selected, the approach does not use the inclusion probabilities 
in the estimation. Under the Bayesian model (2.1) and (2.2) and squared 
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error loss, the Bayes estimator of the finite population mean, Y = J2iLi ^/-^j 
is given by 

(2.3) F B = E[Y\ys] = fnVs + (1 - fn)E[i;' {e)\ys], 

where = n/N and = J2ies ^ the sample mean. 

In order to understand the relationship between the Bayes estimator 
and Ijs for a large sample, we find the mathematical limit of Elijj' {6)\ys] as 
n — >■ oo as a function of y^. Hence, in the following yg is treated as a non- 
random argument of E[ip' {6)\ys]- Consequently, let Tn be a random variable 
with density (as a function of t) 

(2.4) exp[—n(j)~^'ilj{t)]exp[(p~^nyst]^ J exp[—n(j)~^'ilj{T)]exp[(j)~^nysT]dT. 
This leads to the following lemma. 

Lemma 2.1. In addition to (2.1) and (2.2), assume the following regu- 
larity conditions: 

(R.l) The functions V''(-) exp[-(r/2){/i(-) -,9}^]/i'(-) and exp[-(r/2){/i(-) - 
/3}^]/i'(-) are bounded. 

(R.2) The sequence of random variables {Tn} converges in probability 
with respect to (2.4) to T as n — > oo. 

Then E[ip' {6)\ys] converges to C ■.= ijj'{T). 

Remark 2.1. In general, T is a function of ys- This reveals the large- 
sample behavior of Elip' {9)\ys] as a function of ijs- In the Appendix, we 
verify that assumptions (R.l) and (R.2) are satisfied for three well-known 
distributions involving Gaussian, Bernoulli and Poisson distributions, where 
C = ys. 

The following theorem is a simple consequence of Lemma 2.1 and (2.3). 

Theorem 2.1. In addition to the conditions of Lemma 2.1, suppose that 
fn^ f for some < / < 1 . Then 

f^^fy, + {l-f)C 

as Qo. 

In the Appendix, we note that the convergence proof of Theorem 2.1 does 
not use any assumptions regarding the sampling design. In general, for large 
n, is not design-consistent, except possibly for a self- weighting sampling 
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design. Let be any design-consistent estimator of Y. For example, we 
can choose yy^ to be the well-known Hansen-Hurwitz or Horvitz-Thompson 
estimator. Then, using Theorem 2.1, we can obtain the following design- 
consistent estimator of Y based on the Bayes estimator: 

f = Y^-{fy^ + (l-f)C-y^}. 

For a self- weighing design, yw = ys, and so the Bayes estimator is indeed 
design-consistent for the three examples given in the Appendix. 

3. Hierarchical Bayes estimator and its limiting behavior. In this section 
we extend the results of Section 2 and consider two important cases: case (i) 
(f> is known, but /5 and r are unknown, and case (ii) all of the hyperparame- 
ters are unknown. To treat these two cases, we need further information on 
the finite population that allows estimation of unknown hyperparameters 
through an appropriate hierarchical Bayes method. To this end, consider a 
finite population divided into m strata. Let Yij denote the value of the jth. 
observation in the ith. stratum (i = 1, . . . , m; j = 1, . . . , Ni). We consider the 
estimation of a particular stratum mean. Without loss of generality, we con- 
sider estimation of the mth stratum mean given by Ym = N~^J2f=i^mj, 
where Nm is the known population size for the mth stratum. We assume that 
the sample ys consists of rii observations from the ith stratum, 1 <i <m. 

In addition to (2.1) and (2.2), for each Yij (i = 1, . . . ,m;j = 1, . . . ,Ni), 
we assume that f3 and r are independent, with (3 ~ U{—oo, oo), an improper 
uniform distribution over the real line, and r ~ G{a, b), a gamma distribution 
with density proportional to r^''/^)"^ exp[— (a/2)r]. 

Writing 9 = [^i, . . . , 9m]' , we obtain the conditional density 

m rii 

7r(0,/3, r, 4)\ys) oc H H ^w[(l)~^{yijGi - i^iOi)} + p{yij,4>)] 
i=ij=i 

m 

(3.1) X l[{exp[-ir/2){hie,) - f3}']h' {e^y^'} 

i=l 

X exp[-ar/2]r(^/2)~^ 
First, consider case (i). In this case, we show that (see the Appendix) 

E[i;'{em)\ys] 

(3.2) 

= Et„„ [i^'iTnJh'iT^Jg{TnJ]/ET„^ W {T^JgiT^J], 

where the expectation in (3.2) is taken with respect to a random variable 
with density (as a function of t) 

(3.3) exp[(/)"^nm{ymt - ?/'(*)}] / / exp[(/)"^n,„{ymr - '0(t)}] dr, 
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and where g is a function defined in (A. 7). 

This leads to the following lemma, which, in turn, proves the convergence 
of the hierarchical Bayes estimator to ym. 

Lemma 3.1. Assume the following regularity conditions. 

(R.3) The functions ip' {■)h' {■)g{-) and h'{-)g[-) are hounded. 
(R.4) The sequence of random variables {Tn^ } converges in probability 
[with respect to (3.3)] to T as Um — ^ oo. 

Then E[ip' {9m)\ys] converges to Cm = i^'(T). 

Remark 3.1. As before, we can check that the regularity conditions 
(R.3) and (R.4) are satisfied for all three examples in the Appendix and 
that in each case Cm = 27m- Verification of (R.4) is similar to that of (R.2). 
Verification of (R.3) may be facilitated by noting that 

/■m— 1 f rii ^ 

n h'{e,) n eMr'{y^A - m)}] \ 
i=i I j=i J 

m -i-(b+m-l)/2 m-1 

+ Y.h\ei) n 

1=1 J i=l 

For example, in the case of normal distribution, 

„m — 1 rii m—1 
j=l j=l i=l 

The following theorem is an immediate consequence of Lemma 3.1. The 
theorem suggests a simple adjustment to the hierarchical Bayes estimator 
to achieve the design-consistency property. 

Theorem 3.1. In addition to the regularity conditions of Lemma 3.1, 

assume that := Uml^m —* fm for some < /m < 1- Tet be the 

hierarchical Bayes estimator ofYm- Then 

Yjyi, ^ fmym ~l" (1 fm)Cm 

as Um. — > OO. 

Let ymw be any design-consistent estimator of Y m- Then Y m = ^™ — 
{fmVm + (1 — fm)Cm — Vmw} IS desigu-cousistent, based on the hierarchical 
Bayes estimator. 

Now consider case (ii). Suppose (f), jS and r are mutually independent. 
Furthermore, let u := 1/0 ~ [/ (0, oo), /3 ~ U{ — oo, oo) and r ~ G{a, b). 
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In this case, we wih consider the normal distribution only, since (j) is known 
for the binomial and Poisson examples. Since h{9) = 9, integrating (3.1) with 
respect to /?, r and v (in that order), we obtain 

(3.4) 

>, -(nT/2+l) |- rn 1 -(6+m-l)/2 



OC 



m Ui 



<i=lj=l 



a + Y,{Oi-e}' 

1=1 



where ut = Si^i the Appendix, we show that 

(3.5) E[i^'ie^)\ys]^ym. 

4. Uncertainty measures. A conventional measure of uncertainty of the 
hierarchical Bayes estimator is simply the posterior variance V{Ym\ys). 
Noting that V{Ym\ys) = E[{Y^ — Ym)'^\ys], we may define a measure of 
uncertainty (MU) of Ym as 

MC/W = E[{f^ - Ymf\ys] = V{YM + [^m, - F™]^. 

Thus, in order to achieve design-consistency, we increase this measure by 

[Yfn — y™]^. However, this apparent increase may be misleading, since this 
will only happen if the assumed hierarchical model holds for all units of the 
finite population, an assumption hard to justify for the unobserved units of 
the finite population based on the observed units in the sample. 

We now propose an alternative measure of uncertainty following the work 
of Prasad and Rao [20] , who considered a design-consistent pseudo-EBLUP 
for a nested error regression model. Following their approach, we define the 

mean squared error of Ym as 

MSEifm)=E[fra-Vm?, 

where 77^ = E{Ym\Om), the expectation being taken over the first two levels 
of the hierarchical model. Unlike the previous approach, this approach does 
not require extensive model assumptions regarding the unobserved units of 
the finite population, except for the mild assumption of the existence of a 
random effect 6m- This is certainly an advantage of this approach over the 
previous approach. 

Let f]^ = fj^ {(j)) be the Bayes estimator of r/^ . Note that 

MSE{fm)=him{(t>) + h2m{(t>), 

where him{4>) = E[fi^ - 'qmf and h2m{4>) = E[Ym - fjm?- It is possible to 
write down an explicit expression for him{4>)-, although it may not be in 
closed form. 
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Let (p be any commonly used consistent estimator of (j). For example, in a 
mixed linear normal model, we can consider the residual maximum likelihood 
estimators (REML) for the variance components and weighted least squares 
with estimated variance components for the regression coefficients. We can 

estimate by h^mi^) and /i2-m(</') by /i2m((^) = [Ym - V^?, where 

fj^ = fi^{4)), an empirical Bayes estimator of rjm- We propose the following 
as the second niGcisure of uncGrtciinty of ^m,* 

M[/(^) = /ii„(<^) + /i2™(<^). 

(2) 

Note that MUm does not incorporate the variability due to the estimation 
of (j). On the other hand, MUn^ incorporates all sources of variability, but 
may be sensitive to the specification of the prior distribution on (p. Thus, as 
a compromise, we propose the following measure of uncertainty: 

(3) ^ MU'^^ + MU^I^ 

^ m 2 

5. A Monte Carlo simulation. In this section, we conduct a Monte Carlo 
simulation study to compare the performances of MU^^\ MU^"^^ and MU^'^\ 
proposed in Section 4. In particular, we study the performances of these 
measures in estimating the design-based mean squared error defined as 

MSE diY m) = EdiYm — Ymf', where Ed denotes an expectation with re- 
spect to the sampling design. 

We consider 100 finite populations, each of size 60. The main variable 
Y is generated for each unit of the finite populations using a nested error 
model Yij = ^ + Vi + eij, where is the fixed effect, vt ^ N{0, a^) are the 
random effects and eij ~ iV(0, o"g), with the {viYs and the pure errors {ejj j's 
assumed to be independent, 1 <i < 100, 1 < i < iVj = 60. We set /i = 50 and 
ay = 1 and consider two different values of a^, namely (7e = l and 2. 

We draw a sample of size n from each finite population using a proba- 
bility proportional to size with replacement (PPSWR) sampling design. We 
consider three different choices of n, namely 10, 20 and 30. The size mea- 
sure is generated for each unit of the finite populations using an exponential 
distribution with mean 1. 

Flat priors on the hyperparameters are used to obtain the posterior mean 
and the posterior variance needed to compute MU^^^ . We use PROC MIXED 
in SAS to generate 1050 observations from the posterior distributions, but 
only the last 1000 observations are retained for approximating the poste- 
rior means and variances. Essentially, PROC MIXED uses a Markov chain 
Monte Carlo (MCMC) technique. PROC IML is used to obtain the required 
posterior means and posterior variances. For the second measure 
is estimated by the residual maximum likelihood (REML) method. 
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Table 1 

ARB and ARRMSE of the three different measures of uncertainty 



n 




ARB 






ARRMSE 












arrmse'^^^ 


ARRMSE^'^^ 


ARRMSE^^^ 


10 
20 
30 


1 

0.001 
0.011 
0.017 

2 


-0.051 
-0.004 
0.010 


-0.025 
0.004 
0.013 


1.132 
1.340 
1.515 


0.523 
0.673 
0.761 


0.743 
0.891 
1.003 


Oe = 

10 
20 
30 


0.052 
0.045 
0.044 


-0.109 
-0.060 
-0.022 


-0.029 
-0.007 
0.011 


1.266 
1.471 
1.566 


0.726 
0.821 
0.883 


0.935 
1.051 
1.101 



The relative bias (RB) and relative root mean square error (RRMSE) 
for the i-population (1 < i < 100) using the kth measure of uncertainty 
(1 < A; < 3) are defined as 

^^(fc) ^ j^,[M[/f^]-M5E,(F,) 
MSEdifi) 

and 



RRMSE 



MSEdiYi) 



respectively. Table 1 reports the average RB {ARB^''^;l<k<3} and aver- 
age RRMSE {ARRMSE^'''^;! <k<3} for different combinations of (n,cre), 
the average being taken over all 100 finite populations. In terms of ARB, it 
appears that the measure MU^^^ has a slight tendency to overestimate the 
design-based MSE, whereas MU^"^^ has a slight tendency to underestimate. 
This is probably due to the fact that MU^^^ attempts to incorporate all 
sources of variability, while MU^"^^ does not incorporate the variability in 
estimating </). It is interesting to note that for all the measures, ARRMSE 
increases with the increase of n. This behavior can be explained by the pat- 
tern of the inclusion probabilities induced by our sampling design; see [15] 
for details. In terms of the ARRMSE, Mf/^^^ is better than both MC/(^) and 
Mf7(3). Overah, the measure is a compromise between MU^^^ and 



6. Concluding remarks. In this paper, we examine a useful asymptotic 
behavior of the hierarchical Bayes estimator of a finite population mean. This 
leads to a simple method for constructing a design-consistent model-based 



10 



p. LAHIRI AND K. MUKHERJEE 



estimator of a finite population mean. The method is general, in that it can 
be easily applied to both normal and nonnormal cases and is applicable to 
any complex weighting scheme. We have also addressed the important issue 
of measuring uncertainty of our proposed estimator. The simulation study 
reveals that our second measure suffers from a slight downward design-based 
bias. In the future, the Taylor series or a parametric bootstrap method as in 
[14] may be considered in an effort to reduce the bias. It is conceivable that 
our method extends beyond the stratified sampling design, for example, the 
multi-stage sampling design considered in [7] and [17], but this needs further 
research. 

APPENDIX 

Proof of Lemma 2.1. First we compute EIiIj' {9)\ys]- Using (2.1), (2.2) 
and the fact that if a random variable h{9) has density d{-), then 6 has 
density d{h{9)}h' (9) , we obtain the conditional density 

7T{9\ys) a exp[-ir/2){h{9) - P}^]h'i9) 

(A.l) 

X exp[-n(/) V'(^)] X exp[(/) nys9]. 

This yields 

J 'ip'{t)exp[-{r/2){h{t) - I3}^]h'{t)exp[-n(j)-^i;{t)]exp[(j)-^nyst]dt 

J exp[-{r/2){h{t) - (3}'^]h'{t)exp[-n4)-^'4!{t)]exp[(p-^nyst]dty 

which is a function of n and ys- 
Therefore, 

E[i;'{9)\ys] = {i?Tj^'(r„) exp[-(r/2){/i(r„) - P}']h'{Tn)]} 

(A.2) 

- {£;Tjexp[-(r/2){/i(r„) - /3}2]/j'(r„)]}, 

where the expectation in (A.2) is taken with respect to a random variable T„ 
with density (2.4). The proof now follows by using the bounded convergence 
theorem on both the numerator and the denominator of (A.2). 

Verification of the conditions of Lemma 2.1. When (j) is known in (2.1), 
we verify the conditions of Lemma 2.1 for the Gaussian, Bernoulli and Pois- 
son distributions. Verification of (R.l) is trivial. 
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Example 1. Suppose that conditional on ^, y is normal with mean 6 
and variance cr^. Then equation (2.4) becomes 

exp[— nt^/(2cr^)] exp[tnys/c7^]y/ J exp[— nr^/(2(7^)] exp[rnys/(7^] dr, 

a normal density with mean and variance a"^ /n. Hence, T = ys and C = tJs- 

Example 2. Suppose that conditional on 9, Y is Bernoulli with success 
probability 7 = e^/ (1 + e^). The numerator of (2.4) becomes 

e*S-"/(l + e*)" = [eV(l + e*)]"^- [1 - eV(l + e*)]"-"?^^ . 

Hence, Vn ■= e^"/(l + 6^") has a Beta distribution converging to y^. There- 
fore, C = V'(r) = y,. 

Example 3. Suppose that conditional on 9, Y is Poisson with suc- 
cess rate A = e^. The numerator of (2.4) becomes exp[— n(e* — tyg)]. Hence, 



Vn '■= C"^" has the Gamma density e 



nv^nys-l 



I{v > 0)/T{nys). Since £;(K 



ys and Var(y„) — > 0, Vn converges to ys- Therefore, C = tp'{T) = y^. 

Verification of (3.2). First, integrating (3.1) with respect to /?, we 
obtain 

m rii 



(A.3) 



i=ij=i 

m 

X ]l{h'i9i)}r"'^^ X exp[-ar/2]7 



i=l 



X / exp 



-{r/2)Y^{h{e,)-(3y 



.(6/2)-l 



dp. 



i=l 



Then writing ET=i{H9t) " = ET=i{HGi) - h}^+m{(3 - hf, the last in- 
tegral in (A.3) is proportional to r"-*^/^ exp[— (r/2) X)E^i{^(^«) ~ ^}^]- Hence, 



m rii 



(A.4) 



TT{e,r\y,) oc n n ^MrHviA - m)}] X Y[{h'{9i)} 

i=ij=i 



i=l 



X^{(m+fe-l)/2)-l}g^p^ 



(r/2) 



a + Y.{h{9.^-h}' 



i=l 



Integrating (A.4) with respect to r on (0, 00) and using the formula for the 
gamma integral, 

m rii rn 

TT{e\ys) « n n ^Mr^yrA - m)}] x H ^'(^^) 

i=lj=l i=l 
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a+j2{hm-hy 



i=l 



-{6+m-l)/2 



Integrating (A.5) with respect to ^i, . . . , Ojn~i, we obtain 

(A.6) TT{e^ni\ys) oc Y[ exp[(/>~Hy„j6'^ - ^(6*™)}] x /i'(6l„) x g{9„ 



where 



(A.7) 



m — 1 



1 -(fe+m-l)/2 m_i 



n 



1=1 



Note that the function g does not involve (which will be allowed to g 
to infinity). Then, from (A.6), 



E[ip'{9m)\ys]= J ^'{t)T,{t\ys)dt 

= Et„^ W{TnJh'{TnJg{TnJ]/ET^^ [h' {T^JgiT^J]. 
Verification of (3.5). Note that 
(A.8) E[tlj'{9^)\ys] = I evn,Ms)d9/ J VnMs)d9, 

where 



. -(nT/2+1) ( m 

•' U=i j=i J I i=i 

( nra -\ -(nm/2 + 1) 



-{b+m-l)/2 rn-l 
i=l 



m mj - yni) \ 

3=1 ) 

- ymf + ^{y-mj - yrriY 

/ ( m rii 



(n„/2+l) 
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m ^ (fe+m-l)/2 



X < a-' • 



4=1 




-(nm/2+l) 



= |"-m(6'm - ym)^ + - ym)^ j X l{0,n,ys), Say. 

Write KMs) = fc(2/.){n^(^m-ym)' + E"=i(ymi where 
A:(ys) is such that hn^ is a density. Then, from (A. 8), 

(A.9) E[^'{e^)\ys] = I ei{e,ys)KMs)de/ 1 i{e,ys)KMs)de. 



Next, we show that is the p.d.f. of a random variable that converges in 
probabihty to ym as rim tends to infinity. Then the result that the hierar- 
chical Bayes estimator converges to ym. follows by noting the boundedness 
of the function 6 =^ 6l{6, yg). 

For simplicity, write n for Um- Note that /i„ is a density of the form 
c/{n{x — /i)^ + d}"/^"*"^. Clearly, the mean is fi. To prove convergence in 
probability, we next show that for any ei > 0, 

{nx^ + dx/ j [nx^ + dx ^ 0. 

Substituting x = {d/nY^'^y, the above becomes equals /^/2g(y^ + 
/o°°(y^ + where e := ei/d^^"^ . The numerator is bounded above 

by ]^/2^{v}''^ey + l)-("/2+i) dy, which is 0{l/ {n^/'^{n^''^e)'^}). The denom- 
inator is bounded below by /(2e^)~"d7/ = 0{l/n2^). Hence, (A.IO) follows. 
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